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THE METHOD OF CHARACTERISTICS FOR THE DETERMINATION OF 
SUPERSONIC FLOW OVER BODIES OF REVOLUTION 
AT SMALL ANGLES OF ATTACK 1 

By Antonio Feebi 


SUMMARY 

The method of characteristics has been applied for the deter- 
mination of the supersonic-flow properties around bodies of 
resolution at a small angle of attach. The system developed 
considers the effect of the variation of entropy due to the curved 
shock and determines a flow that exactly satisfies the boundary 
conditions in the limits of the simplifications assumed. Two 
practical methods for numerical calculations are given. 

INTRODUCTION 

For the determination of aerodynamic properties of bodies 
of revolution at supersonic speeds, two methods have been 
used: a method that uses the small-disturbances theory and 
a method that uses the characteristics theory. Both methods 
are successful in the determination of the flow properties for 
bodies at zero angle of attack, but the precision of the small- 
disturbances theory decreases when a body of revolution at 
an angle of attack is considered. 

For bodies of revolution having supersonic flow every- 
where, the theory of characteristics can also be used at an 
angle of attack. 

The method of characteristics for the determination of the 
flow field around bodies of revolution at an angle of attack 
was first used by Ferrari (reference 1) in 1936. Ferrari 
considers the flow as potential flow and develops a method 
for the analysis of the flow field around a body that in the 
approximation of potential flow appears to be general and 
can he applied to bodies of any shape and with any angle of 
attack. In the determination of the flow properties along 
the first characteristic surface from which the analysis starts, 
however, Ferrari analyzes the flow around a cone of revolu- 
tion, and in this part of the analysis only small values of 
angle of attack are considered. 

Sauer in 1942 (reference 2) considers the same problems 
and shows that, for small values of angle of attack, the 
analysis of the flow field around a body of revolution can be 
made by applying the characteristic's method .only in one 
meridian plane; and, therefore, Sauer uses characteristic 
lines in place of the characteristic surfaces considered by 
Ferrari. Sauer, in the development of his system, is inter- 
ested essentially in the analysis of the flow around circular 
cones; and when the method is applied to bodies of revolution 
of shapes different from cones, the boundary conditions are 


no longer satisfied. The flow obtained from the solution 
used, also at small angles of attack, wets a body that is 
not a body of revolution. The body can he obtained from 
the body of revolution considered initially by curving its 
axis of symmetry. Sauer also assumes that the flow is 
potential flow. With this assumption, the flow must be 
considered as potential flow for the case of the body at zero 
angle of attack also; therefore, all the effects of entropy 
gradients are neglected. 

The flow field around circular cones at small angles of 
attack has been analyzed in a more exact form by Stone and 
Ferri. (See references 3 to 5.) In his analysis, Stone 
considers the flow as rotational flow and, therefore, takes 
into account the effect of entropy gradients on the velocity 
distribution. This effect exists only when the cone has an 
angle of attack and, at low Mach numbers, is small but of the 
same order as the effect of other parameters that are con- 
sidered in the analysis. Ferri considers correctly the 
entropy distribution at the surface of the cone and found the 
existence of a vortical layer at the surface of the cone across 
which the pressure does not change although density and 
velocity change. 

Here, the method of characteristics is extended to the 
analysis of the flow field around a body of revolution at small 
angles of attack for the case of rotational flow. The effect 
of entropy gradients about bodies of revolution even at small 
angles of attack can be important because the entropy 
gradients that exist in the stream for small angles of attack 
are due to the variation of curvature of the shock existing 
at zero angle of attack also, together with the fact that the 
shock surface does not have axial symmetry with respect 
to the direction of the undisturbed velocity. 

The method presented permits the determination of a 
flow that in the assumption of small angles of attack exactly 
satisfies the boundary conditions and, therefore, wets the 
body of revolution considered. This method is given in a 
form that permits its application to practical problems and 
requires either numerical or numerical and graphical calcu- 
lations of the same type as the calculations used for the 
analysis of the flow around bodies at zero angle of attack. 
The method can be applied to cases in which the entropy 
variations can be neglected or are zero. In these cases the 
terms that contain the entropy variations become zero. 


• Supersedes NACA TN 1809, “The Method ol Characteristics tor the Determination of Supersonic Flow over Bodies of Bevolntion at Small Angles of Attack" by Antonio Ferri, 1849. 
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SYMBOLS 

V , 9 

cylindrical coordinates (fig. 1 ) 

r, 6 

polar coordinates (fig. 8) 

V 

local velocity (function of x, y, 6) 

U, t), w 

velocity components in cylindrical coordi- 
nates ( u along x-axis, v along y-axis, and 
w normal to meridian plane) 

v„ V H , V) 

-velocity components in polar coordinates 
(», along r, v n normal to r in meridian 
plane, and w normal to meridian plane) 

V, 

limiting velocity corresponding to adiabatic 
expansion to zero pressure 

P 

pressure 

p 

density 

y 

ratio of specific heats 

a 

speed of sound (a 2 =y 

a 

angle of attack of body 

0 

Mach angle ^sin 0—^^ * 

<t> 

angle between velocity V and x-axis 

n 

angle between the axis of the cone tangent 
to the shock and the axis of the body 

c r 

angle at the apex of the cone tangent to the 
shock 


tangents to the characteristic surfaces in 
the meridian plane 0 — Constant 

Vs 

velocity component normal to the shock 

- 

surface 

Dt 

velocity component along the generatrix of 
the cone tangent to the shock 

w 

velocity component tangent to the cross 
section of the cone tangent to the shock 

n 

angle between the tangent to the shock 
and the axis of the body 

AS 

entropy variation for unit mass 

n 

normal to the streamline in the plane 
d — Constant 

N 

normal to the surface of the shock 

H, L, K, Z 

coefficients defined by equations (24) 

.4 

coefficient defined by equation (41) 

Q 

coefficient defined by equation (45) 

U\, Dj 

coefficients defined by equations (55) 

•4ij -4ai P u Pi, 

T coefficients defined by equations (60) 

R 

Subscripts: 

radius of the hodograph diagram 

0 

free-stream flow quantities 

1 

flow quantities for the condition of zero 
angle of attack 

2 

flow quantities related to the effect of 
angle of attack as defined in equations 
(5) and (11) 



Figubs 1.— Cylindrical coordinates and corresponding velocity components. 

The prime (') represents quantities in front of the shock 
and the double prime (") represents quantities behind the 
shock. 

EQUATION OF MOTION FOR FLOW AROUND A BODY OF 
REVOLUTION AT A SMALL ANGLE OF ATTACK 

Consider a cylindrical coordinate system in which the 
x-axis is coincident with the axis of the body of revolution, 
the y-axis is normal to the x-axis in any meridian plane, and 
the position of every meridian plane is defined by the angle 
8 measured with respect to the meridian plane that contains 
the direction of the undisturbed velocity (fig. 1).' 

Euler’s equations of motion for steady flow in cylindrical 
coordinates are 


1 bp bu bu . bu 
bx~l)x U +by V ~^yb9 


P 

1 bp bv 


bv 


■ . bv 

P by bx U ^ by v ybd 

1 bp bw bw bw 
P ybd'~bx U 'by U ~^yd£ 


The continuity equation in cylindrical coordinates can be 
expressed in the form 


■w 

(la) 

w 1 
to— — 

y 

(lb) 

. vw 

r«H 

? y 

(lc) 


b(pu) . b(pvy) d(w/>) 
5x y by. ^ ybd ~ 


=0 


( 2 ) 


whereas the law of conservation of energy can be written in 
the form 


/I dp 

V 2>p\ , 

( & u . 

bo. 

bw\ 

(3a) 

\P dx‘ 
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(“ 81 + " 
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/ bu , 
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If the density is eliminated from equation. (2) by means of 
equations (1) and (3) and the quantity a is introduced as 
defined by 



Between rotation of flow and entropy the following relation 
exists: 

curl VXV=^ grad S (7) 


or for small angles of attack 


the following equation can be obtained: 



In this analysis only small angles of attack will be con- 
sidered, and, therefore, only the first-order effect of the 
angle of attack will be determined; whereas, the quantities 
of the same or higher order than the square of the angle of 
attack will be neglected. In this approximation the velocity 
components of the flow around the body can be expressed 
in the following form (references 1 to 4) : 


u=Ui+au 3 cos 6 

(5a) 

v = v 2 J rav 3 cos 8 

(5b) 

w=aw 1 sin 8 

(5c) 


where u, v, and w are functions of- the three coordinates x, y, 
and 0; whereas u u v u u 2 , % and w 2 are functions only of the 
coordinates x and y of any meridian plane. The .quantity a 
is the angle of attack of the body, the quantities with sub- 
script 1 are the quantities existing at the position (x, y) for 
the body considered at zero angle of attack at the same Mach 
number, and the quantities with subscript 2 are functions 
that take into account the effect of angle of attack. 

It will be shown in the following considerations that the 
form assumed in equations (5) for the velocity components 
permits the boundary conditions to be satisfied in the simpli- 
fications assumed. For small angles of attack equation (4) 
becomes 


du 

bx 


/ u*\ bo/ v*\ v /bu , baXur 

V 1 a*r?yyV a 2 ) + y \by^bx) a 2 ' 


bw 


ybe 


= 0 (6) 


Equation (6) is similar in form to the corresponding 
■ equation for the case of the body at zero angle of attack and 

differs only in the term In order to analyze the differences 

between this expression and the expression for the axial 
symmetrical case and in order to obtain another relation that 
defines the quantity w, the relation between rotation of the 
flow and entropy gradient will now be introduced. 


bS a 2 / bv bu\ 
bx yR~ V \bx by) 

bS a 2 /bv bu\ 
by yR~~ U \bx by) 


(8a) 

(8b) 


bS a 2 

(bw bv \ 

l ] ui 

/ bu bio\ 

y b8 yR V 

K by y be) 

1 ' Ur \ 

y 

\jybd bx) 


(8c) 


If n is the nor mal, in the meridian plane 0= Constant, to 
the local tangent to the streamline, then 

bS bS v . bS u_ 

bn bx 'V by V 



whereas from equation (8c), when equations (5) are used, 
it follows that 


1 a* bS 
y y R bfi 


r„ 

v 3 \ 


/Ui 

bw 2 \~ 


yf 

-f- \‘ U, 1 

y 


bx/_ 


( 10 ) 


Therefore, the entropy can be expressed in the form - 


AS=AS l J r a.AS 2 COS 0 (11) 


where AS t and AS 2 are functions only of x and y. Equation 
(11) is valid outside of the vortical layer existing at the 
surface of the body (see reference 5); whereas the entropy 
at the surface of the body is constant and equal to the- 
value existing at the conical tip. Because no pressure 
gradient exists across the vortical layer, the pressure will 
be determined outside of the layer and equation (11) will 
be used in all the flow field. The presence of the layer must 
be considered for the velocity and density distribution at 
the surface of the body. At the surface of the body the 
entropy is constant and equal to ASi — a ASj. (The meridian 
plane 0=0 is defined as in figure 1.) From equation (11) 
It results 


v 


bw 2 . bw i 

■W +u -sr 


VW 3 -\-VVi+'WU 2 

y 


yyR 


AS, 


( 12 ) 


Equation (6) can be written in the following form: 

bu / u*\ . bv / v 2 \ v 2uv bv . uv bS , bw 

bx \ ffiV"'dyV ay 'y a 2 bx VyRbn ybd 

( 13 ) 



1042 


REPORT 1044 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Equation (13) together with equation (12) defines the law 
of motion of the flow around the body at small angles of 1 
attack. These equations will be used as a basis for the cal- 
culation of the flow field by the method of characteristics 
to be treated in a subsequent section. 

CONDITIONS AT THE SHOCK FRONT 

Equations (5) and (11) represent a stream that wets a 
body of revolution at a small angle of attack. In order to 
satisfy the boundary conditions at the surface of the body, 
the functions u 2 , v 2 , and w 2 must be properly selected. 
Equations (6) and (11) must, however, satisfy the boundary 
conditions at the shock surface also in order to be a solution 
of the problem. It is necessary, therefore, to show that a 
shock surface can exist across which the undisturbed stream 
inclined at a with respect to the axis of the body is trans- 
formed into a flow represented by equations (5) and (11). 

In order to show that the shock boundary conditions can 
be satisfied, the following procedure will be employed. A 
shock surface distorted in a maimer to be described is 
assumed. Then, the free-stream velocity ahead of the shock 
will be resolved into three components: normal to the 

shock in the plane 0= Constant, v T tangent to the shock in 
the plane 0=Constant, and w perpendicular to the plane 
0= Constant. Similarly, the flow behind the shock will be 
resolved into three components. In addition, each com- 
ponent of the flow behind the shock will be divided into two 
terms: one term for zero angle of attack and one term for 
the difference due to the angle of attack (for example, 
«=u 1 +att i cos 0)'. Then, the conditions of equilibrium at 
the shock will be imposed, and it will be shown that the 
terms U 2 , v 2 , and w 2 at the shock are independent of 0 when 
the angle of attack is small as initially assumed; hence, the 
distorted shock is consistent with the field of flow behind it. 
Such a shock surface can be obtained by deforming the shock 
surface produced by the body when the angle of attack a 
is zero in the following way (fig. 2): 

When a=0, the shock surface is a surface of revolution 
in axis with the body; therefore, if OP, OP' is the curve 
intersection of the shock with the meridian plane 0=0, 
then for a=0 the tangent AQ at any point Q of the curve 
OP is the generatrix of a circular cone having the vertex at 
a point A of the axis and tangent along the circle QQ' to the 
shock surface. The shock surface, therefore, can be consid- 
ered as a surface envelope of circular cones having the axis 
coincident with the axis of the body but having variable 



cone angle and variable position of the apex A along the 
axis AB of the body. For the case of a ^ 0 the shock surface 
is not a surface of revolution but can still be considered, 
for small angles of attack, as the envelope of the same 
circular cones considered for the case a=0. These cones 
have the same apexes and the same cone angles as the cases 
for a=0 but do not have the axis of symmetry AB coincident 
with the axis of the body AB although they are rotated in 
the plane 0=0 with respect to the body axis. The angle 
V, through which each axis of the cones must rotate in the 
plane 0=0 with respect to the axis of the body, is not 
constant but varies for each cone considered. For example, 
the cone AQQ' tangent to the shock surface for a=0, when 
0, must be rotated by an angle i; to the position AQiQ/; 
the axis ABi remains in the plane 0=0. 

The shock surface so generated is consistent with tho flow 
represented by equations (5) and (11), and this can be shown 
in the following way: 

Consider a point P of the shock produced by the body at 
an angle of attack, and consider the cone tangent to the 
shock at the point P (fig. 3). Call o- the semiangle of the 
cone with respect to its axis of symmetry. The axis of tins 
cone is inclined at an angle ij with respect to the axis of the 
body and lies in the plane 0=0. 

The unif orm velocity V 0 ahead of the shock is decomposed 
in the three components: v N ' in the direction PB normal to 
the shock, v r ' in the direction AP along the generatrix of the 
cone, and w' in the direction normal to the plane APB. 
These components are, at small angles of attack, 

Vif'—Va sin o'— F 0 (a— ij) cos a cos 0 (14a) 

V T '=V 0 cos a -\-Voia — )?) sin a cos 0 (14b) 

«/= — Vo(«— ij) sin 0 (14c) 

Strictly, in equations (14) ^ must be written in place of 0; 
but, for small angles of attack in equations (14), the differ- 
ence between and 0 can be neglected. Indeed, 

I 

^= 3 + 0 ' 

where 3 is of the order of a, and 0' differs from 0 by a quantity 
of the order of a. 



Figtjbe 3— The velocity components In front of and behind the shock. 
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The velocity components behind the shock are (fig. 3) 

Vk"=u sin a- — v cos 0-+17 cos 8{v sin <r+-u cos a) (15a) 

v T n =% cos <r+jj sin o-+)j cos d(v cos a— u sin a) (15b) 

w"=w-\-ri(u — v cot <r) sin 8 (15c) 

-where u, v, and w are the velocity components behind the 
shock in cylindrical coordinates in axis with the body at the 
point P considered. 

The velocity components u, v, and w at the point P can be 
■expressed in the form given by equations (5), in which the 
■components and Vi are the quantities obtained at the point 
P for the condition of a= 0 and are direct functions only of 
x and y. The point P, however, is a point of the shock, and its 
coordinates x and y change when the coordinate 6 changes; 
therefore, the velocity components iq and Vi at P also change 
with 0i. In order to separate the part of the components u, 
v, and w dependent on 0 from the part independent of 0, the 
velocity components Ui aDd Vi at P will now be expressed as a 
function of the flow properties at a point P t near P, having a 
constant value of x and y for every value of 0. 

Now, it has been assumed that the angle of the cone a 
tangent to the shock at the point P is equal to the angle of the 
cone tangent to the shock for the condition of zero angle of 
attack at the point Pi (fig. 4). The point P t is obtained on 
the shock by rotating the cone APQ tangent to the shock for 
the condition a=0 through an angle 17 around the axis AN 
normal to the plane 0=0 at the apex A of the cone. Because 
for the condition of zero angle of attack the velocity com- 
ponents Ui and ®i are independent of the coordinate 0, the 
velocity components Ui and Ui at P (x P ,y P ) (fig. 4) are equal 
to the velocity components at Pafe,^) in the plane APiC. 


Therefore, if A N is the distance PjPi, 


- 


(16a) 



(16b) 

where (fig. 4) 

Xp,V 

AN = — — cos 0 
cos 0- 

(17) 


Substituting equations (5), (16), and (17) in equations (15) 
results in 

®jr"=(a 1 sin <r—v 1 cos ff) Pl +a cos 0(ii 2 sin <r— cos a-) Pl + 

17 cos 6(u 1 cos <H-»i sin <r) Pl + 

x Pl y /d«i . di’i \ . 

— - cos 0 ( sin (r — rr. cos x ) (18a) 

cos <r \bN b N Jpj, 

c r /r =(w l cos <r+Vi sin a) Pl +a cos d(u 2 cos tr+» 2 sin cr) Pl — 

17 cos 0(u 1 sin a — a: cos <r) Pl + 

X P L V a/'bUi . diJi . \ 0, S 

— - — cos 0 ( ■^rt cos & sin <r ) * (18b) 

cos <r \bN bN Jp x 

vf=an i/ 3pi sin 0 + 17(0*1 — v 1 cot <r) Pl sin 0 (18c) 



Figure 4 u — Th& velocity components behind the shock foe as^O as a function of the Telocity 
components behind the shock foe «~0. 


For the condition of zero angle of attack at Pi 

Vtr "=Ui sin <r—Vi cos 0- (19a) 

v t "=Ul cos cr+ Wi sin <r (19b) 

and for .the condition of the equilibrium at the shock at zero 
angle of attack 

V 0 cos (T—v t" (20a) 

W" To sin £ r=^j(F, 2 - i7r 1 " 2 ) (20b) 

At the point P for the case of a small angle of attack, 

(i>/^=^(F { *-a/*) P (21a) 

v Tp "=v Tp f (21b) 

w"—w' (21c) 


If equations (19) are used, equations (18a) and (18b) can 
be written in the form 

(s/)p=(5q) Pl +“C 0 S 0 (%/) Pl +lI COS 0 (v T ") Pl + 
x p ,v 

— COS 0 ( ) 

cos <r \ CUV / 

(27r'0i >= (®r 1 / 0 Pl +a COS 0 (vt 2 ") Pi —V COS 0 0>jr 1 ")p 1 + 
cos 0- \ bN / Pl 

Therefore, from equations (14), (20), and (21), 

— (a— q)Vo COS <r {Vtr") Pi fr 


( x P i) bvir‘ 

1**1"+^ ~w 




2(y— 1) 
7+1 


(a— t]) Vo cos <r sin 0- 


(a— i 2 )F 0 sin < T=(vT a ")p 1 «— 17(0*/%,+ 

-I 


x Pl v / bv tI 
sin<r\ bN 


(22a) 


(22b) 
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or 




n gp, l ( ^2lL 
a sin cr V 0 \ bN 


bv r ''\ 
bN ) Pl 


(23 a) 




sin 0 -+ 


(l-l) 


v ■ 

atn 

<r V XPl 


V") 1 

V «) 

\ v n ) Px 

Sill 

CL 

a COS a 

IV 

l £>iv ) P J 


fvh\ _ A v\ i ( —Ui + Vi cot cr \ rj_ 

\yJpT v «/ + V Vo J Pl a 


(23b) 

(23c) 


or for small angles of attack 


tan <j>=- 
u 


and 


sin j3 = 


VP+t? 


The tangent to the line intersection of a characteristic 
surface with the meridian plane 0 = Const ant is 


Xa=tan (<t>+P)—-g — jj VS* — HL 
X 6 =tan (.4>—P)=j^ J rjj ^K*—HL 


(25) 


In equations (23) the coordinate 8 does not appear; there- 
fore, for the shock considered the functions «j, and v>i 
are independent of 6, and equations (5) represent a flow con- 
dition in agreement with the conditions at the shock. 

The ratio r\ja which appears in equations (23) is independ- 
ent of a; therefore, for a given point ij/a remains constant 
in all the range of angle of attack in which the simplifica- 
tions assumed are valid. (Indeed, Ut, v 2 , and io 3 are also 
independent of the angle of attack (equations (5)). The 
values of u a , v lt w a , and ij/a must therefore be determined 
only for one value of the angle of attack. 


METHOD OF CHARACTERISTICS FOR FLOW AROUND A 

BODY OF REVOLUTION AT A SMALL ANGLE OF ATTACK 

Lx this section the method of characteristics is applied to 
equation (13) to establish equations which will permit the 
flow field behind the shock to be calculated by a point-by- 
point process. If the flow is anywhere supersonic, equations 
(12) and (13) permit the determination of the flow around 
a body of revolution at a small angle of attack by using the 
method of characteristics. Equation (13) can be written in 
the following form: 

where 


L 






Z 


uv 1 bS . bw . 

V 7 R bn~^y b8~* 


- 

yj 


(24) 


If <j> is the angle between the velocity V and the z-axis and jS 
is the Mach angle, 


tan <f>= 




where X* is the tangent to a line corresponding to the char- 
acteristic surface of the first family and X» is the tangent to a 
line corresponding to the characteristic surface of the second 
family. The terms X« and X t are solutions of the equations 
(reference 6 ) 

27X 2 -2KX+Z=0 (26) 

Because.^, v, V, and a can be considered to be given by an 
equation of the type of equations (5), <i> and /} can also be 
writ ten Jn the form 


CC(f>2 cos 6 

0=/3i+a£i cos 9 
> 

The characteristic surfaces are not, therefore, surfaces of 
revolution but can be obtained, as was true for the case of the 
shock, as an envelope of circular cones with their apexes at 
the axis of the body and their axis of symmetry hi the plane 
0=0 and inclined with the axis of the body. 

The determination of the u and v components of the 
velocity in any point of the flow can be obtained from equa- 
tion (13) by performing a transformation in order to obtain 
a law of variation along the characteristic lines (reference 6 ). 
Indeed, for every point of any meridian plane (for example, 
of the meridian plane 0 = 0 , or 8—k) two characteristic lines 
can be obtained as the intersection of two characteristic sur- 
faces with the meridian plane. Along these lines the varia- 
tion of the u and v velocity components is determined by the 
equations of characteristics that can be derivod from equa- 
tion (13). Assume that at two points Pi and P 2 (fig. 5) of 
the meridian plane 0=ConBtant (for example, 0=0, or 0 =t) 
the velocity components are known. From equations (25) 
the tangents to the characteristic surfaces in this meridian 
plane can be drawn and the velocity components u and v at 
the point Pj, intersection of the two tangents, can be ob- 
tained in the first approximation. 

The equations of characteristics can be obtained by 
analyzing equation (13) along the characteristic lines given 
by equation (25) in the following way: If du and dv are the 
variations along the characteristic lines, 



du bu . . bu 

dx~bx +( - Xa ’ 


u 
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0-0 



Floras: 5— The analysis of the flow with the characteristics system. 


or (see equation (9)) 




'bx 


a* bS 
yRV bn 


then 


dv 

dx 




dg 

by 


du du dv bv a 2 bS 

dx~bx^ K dx by K yRV bn 


(27a) 

(27b) 


If equations (27a) and (27b) axe substituted in equation (13) 
and equation (20) is used, along the characteristic line of 
the first family defined by 

^=X tt =tan(j3-{-<£) (28a) 

there results 


du . . d v 
dx + ^dx 


a 2 dS 
yRV dn 


(h X “) + (y + ybd)H ~ 0 (28b) 


and., therefore, is known at the points Pi and P 3 . The 

value of the entropy is also known at the points P t and P 2 and, 

therefore, the value of can be determined (reference 6) 

dS_ ASp 2 -AS Pi 

dn' 




(29) 


From equations (28) and (29) the values of u and v can be 
determined in the first approximation for the point P 3 . In 
order to determine the value of w at P* the following Proce- 
dure can be used: 

If 8 is the projection of the streamline in the meridian 
plane considered (fig. 5) and P s is a point near Pi and P 2 , 


dw 2 bw a • , , bw 2 

-^=-^-sm*+ w c° S 0 


bw 2 o 
by V 


bw 2 u_ 
bx V 


or from equation (12) 


bw 2 
bs ~~ 

Now (fig. 5) 

Wt B =W ip + 
W * F = W *P t + 


dw 2 +vv 2 -\-uu 2 sin 2 ft V AS 2 


yV 

y- By 


r cos ft "I 

Lcos (ft + 0)Jp 


r cos ft "I 

Lcos (tf— jS)Je 


(30) 


(31a) 

(31b) 


and 


( dwj \ _ 
dn/ E 


w 1f -w 2b 


(32) 


Therefore, 



and along the characteristic line of the second family defined 
by 

^=X 4 =tan (<£— / 3) (28c) 

there results 


du , , dv a, 1 dS fK \ , /g , 1 

dx +Ka dx yVR dn\H ^rKy^ybeJS 


0 (28d) 


Equations (28b) and (28d) contain the term but at 
small angles of attack 


bw WtO. . w 

— = cos 8 =— cote 

ybe y y 


The values of u, v, w, AS, u u v lt and AS t are known at the 
points Pi and P 2 ; therefore, the values of itj, v 2 , v%, and AS 3 
at the same points can be calculated from equations (5) and 
(11). (The values of u Y , v u and AS t at those points are 
known from the determination of the flow for a—Q.) There- 
fore, from equations (30) to (33) the value of w 2 at P 3 can 
be determined. 

After the velocity components u, v, and w at P 3 have been 
determined in the first approximation, a second approxima- 
tion can be determined by assuming the average values 
between the corresponding values at the points P 2 and P 3 or Pi 
and P 3 for all the coefficients. After the velocity components 
at a point P 3 have been obtained, the velocity components at 
any other point having the same x and y coordinates as 
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P s but a diff erent coordinate 6 can be calculated from 
equations (5). 

For practical calculations, equations (28) can be trans- 
formed in the following form: 



\j=tan (^-j -<f>) 

(34a) 


~Y~ d * taa ^ Lv - ^ ^ y v be ) taiLjS 



sin 2 /3 dS~\ sin (8 , n 

yR in J cos (<£+£) 

s 

(34b) 


X 6 =tan 3) 

(34c) 


Tr+^e 



sin 2 /3d£~] sin fi , 

yR dn\ cos (0 — /3) aX 

(34d) 


1 dw 2 iVi w 2 . sin 2 j3 AS 3 

V ds y V Vy** 1 * yR y 

(34e) 

and 

OC(f>2 COS 6 

(35) 

where 

D 2 • 

^2" COS <f>i ^ 8UL <j> i 

(36) 

and 

V=Vi+aV 3 cos 0 

(37) 

where 

Vi—u 2 cos <j> L +v 2 sin <&. 

(38) 


At the surface of the body outside of the vortical layer the 
calculations are similar to the case of zero a ngle of attack 
because the entropy at the surface of the body is known 
in every meridian plane and the value of 9 is given. Equa- 
tion (34e) gives the variation of w 3 along the body outside of 
the layer; therefore, the value of Wj can be obtained directly 
from another point on the body in the same meri dian plane. 

At the surface of the shock the system of calculations is 
similar to the system for zero angle of attack. In figure 6 
the point P 8 is at the intersection of the tangents to the first 
characteristic surface at Pi and to the shock at P$ in the 
meridian plane 0=Constant. The equations of the shock 
and equation (34b) must be verified at P 8 , which is assumed 
as a point of the shock in the first approximation. 

Tn the plane 0=0, w is zero and the values of V, AS, and 
<f> behind the shock are functions only of the value of £2; and 
for any value of 0, the values of V, AS, and <f> can be 
obtained from the equations of the shock 


cos (Q— g) _ V VfVt 
cos (G — <f>) Vo VofVi 


(39a) 



and 


log. t(~ T ) 7+1 [m„* sin 2 (0-aj>- 
V' ] [M 0 2 sin i (G-£ E ) +:L 2~ ] (39c) 

f 1+ RW C39d > 


If the plane 8—t is considered, the sign of a in equations (39) 
must be reversed. 

From equations (39) the values of V and AS can be deter- 
mined as a function of <j > ; then and as a function of <j> 
can be evaluated. Now, if 4 > P i is the direction of the velocity 
at Pj, the velocity at P 8 will have the direction 

<t>p i —4>p l J rA4> 


Therefore, the velocity at P 8 must correspond to a deviation 

:ai 

dl 


across the shock of 4>p 3 — a and can be expressed as 








where V* p is the velocity behind the shock corresponding 
to the direction In a similar way, 

AS,.=4V 1 +( !? £rX,A* 

Therefore, equation (34b) at the point P r becomes 

K sin 4> , w 2 a cos 8\ sin /3 tan jS“I , 
y yV~ ) cos (P+tiJp, 

K «•> 


1 'yf 1 Afp 2 

tan (<f> — a) _ 2 M* sin 2 (G— «}— 1 


‘] 


tan (G — a) (39b) 


In equation (40), A<f> is the only unknown and, therefore, 
can be determined. From the value of <j> the value of Qp 3 
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and the value of T r p can be determined; a second approxima- 
tion for the position of P s and its value of the velocity can 
be calculated if the corresponding average values between 
P a and Pi are Assumed for <j>, /9, and all the coefficients of 
equation (40). 

The value of w at P can be obtained from equation (23c) 
in which n is given by figure 6 as 

dN . o 
n = — sm 12 

Vp, 

where 77 = 0 — 0 - for 0=0 send 77 = 0-— 0 for 8=r. The value of 
a corresponding to the point P t on the shock for a =0 is 
given by the relation 

Vp* _ Vp* 
sin <r sin O 

and yp i =f(<r^ is the curve that represents the shock for a= 0 

PRACTICAL APPLICATION OF THE CHARACTERISTIC 
SYSTEM 

GRAPHICAL NUMERICAL METHOD 

The analytical part of the characteristic system used for 
determining the flow field about a body of revolution at an 
angle of attack is similar to the system used for a body of 
revolution at zero angle of attack (reference 6 ), but the 
practical numerical application is slightly more involved. 
In equation (34e) the values of Y 2 and A S 2 must be known 
in order to determine the value of w 2 and must be determined 
from equations (37), (36), and (11), where the values of Vi 
and ASi are considered known in the entire flow field and 
given by the determination from the case of zero angle of 
attack. In the practical case, however, the values of Vi, 
<£ 1 , and A&\ have been obtained with the characteristic system 
only in a finite number of points at the intersections of the 
characteristic net, and the characteristic net for the case of 
zero angle of attack is different from the net used for the 
case of a body with a small angle of attack. In order, 
therefore, to obtain the values of V t and A Si at the inter- * 
sections of the characteristic lines for the case with a given 
angle of attack, a lengthy interpolation of the values Vi and 
ASi would be necessary if the two characteristic nets for 
zero angle of attack and for a given angle of attack were 
constructed independently. 

In order to reduce the numerical work to a minimum, the 
two following methods can be used, the first of which is 
practical when a graphical numerical calculation is per- 
formed, whereas the second can be more convenient when 
automatic computing machines are used. 

In both cases the calculations start with the determination 
of the flow at an angle of attack around a cone when the 
body considered is a pointed-nose body of revolution or with 
the deter mina tion of the shock at the lip of the body if the 
body is an open-nose body of revolution. (The tangent to the 
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shock at the lip can be determined with the two-dimensional 
theory.) The flow around a cone at an angle of attack 
has been deter min ed and tabulated in reference 3; whereas 
the flow for zero angle of attack has been tabulated in 
reference 7. A different method for determining the flow 
around a circular cone at an angle of attack is given in the 
appendix. It can be assumed, therefore, that the flow along 
the first characteristic line of the first family at the end of 
the conical region in the plane 0 = Constant (for example, 
8=v) is known (fig. 7). 

For the practical numerical calculations a value of the 
angle of attack must be selected. In order to obtain higher 
precision, it is convenient to select a relatively high value of 
the angle of attack because in this way the differences be- 
tween V and Vi and AS and ASi are large and, therefore, can 
be determined with sufficient precision. 

Usually, when the dete rmin ation of the flow field for the 
case of zero angle of attack is made with a graphical numeri- 
cal process, in order to avoid numerical errors of computa- 
tions, the value of the intensity and direction of the velocity 
are plotted as a function of the position along the character- 
istic lines for both families of characteristic lines. The 
velocity distrib ution and the entropy-variation distribution 
along the characteristic lines and along the surface of the 
body for the case of zero angle of attack can therefore be 
considered known. If the distribution is not given, the 
values of Vi and A Si must be determined as a function of x 
along each characteristic line of a given family (for example, 
of the second family) along the body. 

Then the construction of the characteristic net for the 
selected angle of attack must start by drawing the first 
characteristic line PoPjPj over the design of the characteristic 
net for zero angle of attack (fig. 7). 

From equations (34c), (34d), and (34e) the flow at Pi can be 
determined. From P! and P 2 the point P s can be obtained 
in the first approximation as the intersection of the tangents 
at P 2 and Pi to the characteristic lines. By using equations 
(34b), (34d), and (29), V, <f>, and AS can be obtained in P a 


Characteristic net far a small angle of attach 
Characteristic net for zero angle of attach 
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where 


as for the case of zero angle of attack (reference 6). From 

the variations and along the line PiP s the values 

of V, <P, and AS at the point P 4 can be obtained, where P 4 
is obtained from the intersection of the characteristic line 
PiPj with a characteristic line of the second family in the net 
for zero angle of attack. At the point P 4 , V h AS U and fa are 
known and, therefore, fa, V it AS 3 , and w? can be obtained. 
From the values obtained from the first approximation a 
second approximation can be obtained. From P 4 and P 4 
the point P# can be determined in a similar way, and the 
flow at P 7 can be calculated. By proceeding in a similar 
way, all the flow field can be analyzed. 

NUMERICAL METHOD 

The equation of motion (13) can be transformed by means 
of equations (5) and (11) in a system of equations that 
permits a numerical determination of the quantities V 3 , fa, 
and A St. This system is numerically more involved; 
however, the characteristic net determined for zero angle of 
attack is used. For a small angle of attack, 


where 


7 — 1 

a i 1 

Therefore, 


l 1 r y 1 "I 

-J =— 2 1+-TT (U|U>+t>i*i)«cos 8 

Qi d\ L tfri J 

Yj-= A (41) 


(«,«,+ tqt> a )=^ ViV t ~ 


4=^(1 +Aa cos 8) 

Substituting equations (5) in equations (6) results in the 
following expression if higher-order terms are neglected: 


bUi 

dr 


/. u, l \ rj t>A UiVi /drq . drA . v i 

V <fa*J + dyV <hV ~a? ~\by + *xj + y 

. i dtq f2u i u l -\-‘ifa i A\ , 

L V V 5* \ <h / 

but / ttA dp, / 2a a «i +«i 2 A \ . dv, / _oA_ 
da; \ a,V dy \ <*1* / V «iV 

Ui»i/du s , dPiYI, 40 ^ 
a, 2 \i>y + dx) + ' 

Because the left-hand side of equation (42) must be zero 
for the conditions at zero angle of attack, 

dti s /\ «A , do*/, «A , D, , w 3 tti®i/dtt 2 , dc a \ , 

(43) 


6 


/d4i! 2 

d-»i* N 

\ 2u % +UiA 1 , 

W | 

V dz 4 

dz , 

> fa* 2 1 

Uk + 

4 "du^ 

V U } V 3 — U 2 t?i 


5 dy , 

f a, s 



(44) 


or, from equations (9) and (36), 


0 dVi 2u a +UiA „ bVi 2 V t+ViA faVi bS x 

1 da: a* 1 d y 1 a/ yR dn t 

(45) 

The value of all the coefficients at the points Pi and P a in 

equation (43) can be considered known because and -r~ 

can be considered known from the calculations for the case 
of zero angle of attack. Therefore, equation (43) can be 
considered an equation in which the characteristic linos are 
equal to the characteristic lines for zero angle of attack 

because the coefficients of the partial derivatives 1 ■—} 

and ^ are the same in both cases. Tims, 
da: dy 


(^)=Xi« : = t a n (fa+00 

(jiz)~ tan (fa — fii) 


(46a) 

(46b) 


Equation (43) can be transformed by introducing the 
dS 

entropy gradient > and the equation of motion along each 

characteristic line can be obtained. From equations (9) 
and (11) by means of equations (5) and (37), the following 
relation can be obtained: 

dS d&.dSj . 

a cos 8 

dn on on 


/d Si , d St „\ Si+f a “ cos 8 , 

= -{- d^+- d7 “ C0S V K+V 1 aco- S -8 + 

/dSi . d Si A fa+Uaat cos 8 
(Hy + Tiy “ °° s v U+ cos e 

or, for small angles of attack, 

d S dSi . d& .,r d S lf Tjr T7- \i 
dn = d7h + dnr “ C ° S 0+ L ( ° a Vl ~ Vl Va)+ 


d&l , rr • Tr .”]0!C0i 


COS0 
2 


(47) 
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From equations (36) and (38) 


Ui= P 2 cos fa— Pi sin fa 

(48a) 

v t = V 2 sin fa+faVi cos fa 

and 

»i=Pi sin & 

(48b) 

Ui— Pi cos fa 

Therefore, 


-H (c 2 Pi-5iP s )+^ (UiVt-UiVi) 


. rrs/aSi 1h , bSi «i\ 

- Ux p+ay pj 


Because the term in parentheses on the right-hand side of 
this equation represents the variation of entropy along the 
streamline, which is zero, equation (47) becomes 

as as t , a s 2 

-—=—-±±—1 a cos 0 

UTl UTli OTli 

(49) 

Then, from equations (9), (41), and (49) 


a« 2 a« 2 oi s as 2 /p 2 . ,\ a x * bSi 

ar a y yRVtbn L VP t 1 JyRVi a^ 

(50) 

The equations of motion along the characteristic lines de- 
fined by equations (46) can be obtained by means of trans- 
formations similar to those of equations (27) and are 

diz 2 +Xu, dv 3 +Di dx=0 

(51) 

Xu— ^^=tan (<£i+0i) 

(52) 

dt^+Xia dvi+Di dx=0 

(53) 


In. equations (51) and (53) the coefficients D t and D t 

contain the derivatives and ^ that must be obtained 

ox oy 

from the analysis of the case with zero angle of attack. 
Now, for every point P the variation of Pi along the charac- 
teristic line of the first family for the case of zero angle of 
attack is 



(56a) 


whereas along the characteristic line of the second family 


' dVi\ a P t , , a Pi 

. dx At. a* " t_16 by 


(56b) 


At every point P given by the intersection of two charac- 
teristic lines Xio and Xu in the characteristic net, the values 


fdVA 


and 




\dxj^ Au 


are known, having been obtained 


from the evaluation of the following equations (reference 5): 


dy 

dx 


=Xi a =tan ($ i+/3 


0 


(57a) 


1 / dPA fd(f > A sin fa sin tan 1 

PtVdrA^ t&liPl VdxJ cos Ori+jSi) y 


dS i 1 sin 3 ^ 
dn i yR cos (<£i+j3i) 


(57b) 


dy 

dx 


— X is =tan (fa— Pi) 


(57c) 


+tan p t ( Sin ^ ^ tan * 

Px V dr Au P \dxj y coa(fa— ft) 


dSi 1 sin 3 ^! 
dni yR cos — /3i) 


(57 d) 


Xi» — — tan ($j — 


(54) 


where 


Di= (?£i + »i + q)^L il_ 
\y V 7 ai 2 -Vi s 




ra Si_ 


7 — 1> 

|Mf] 

La?ii 

pA 1 

sin 2 fij 

anj 




(55a) 

rag 2 

v *( 


\ap-| 

La^ 

~viV 

1 sin 1 fay 

f dnj 


(55b) 


Therefore, the values 


and 


a Pi _ Xis / dPj\ Xu / dV A 

ar Xu lia \ dx /Xjj, X IS — X te \ dx /tn 

aPi_ i r/dPA / dPA i 

a y X 16 — X^ Lv dr A« V dr Aid 


(58a) 

(58b) 


can be calculated directly for every point of intersection of 
the characteristic line (equations (57a) and (57c)). 

After substituting the expressions of equations (48) and 
(58) in equations (51) and (53) after some simplifications 
and trigonometric transformations, the following equations 
can be obtained: 

Xu=||=tan (fa+fa) (59a) 


21 S 6 ST — 53 67 
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<iV 2 j. _ a j , , f sin 3 £i 1 bS 3 , , A , 

n ~ tan w+** A '+ 

V 2 B Wi 1 tan fa sin 071 „ 

V, 1 V,* cos'(*+'j8i) A 


X ls =^|=tan (0i— ft)- 


(59b) 

(59c) 


dv 2 

TT 


-tan 


Pldfa — |j 


w 2 1 tan ft, 
Vi y cos (0i 


1 £)5j 
cos (0£— Pi) yR 3?ii 


sin 3 Pi I ooj , . 

<P2-a-2 — y 


>1— Pi) J 


^P,+ 


(59c 


where 

A== i I sm. 

1 COS Pr COS (0i + fti)|_ y 

J_ (dVj\ cos (0i— Pi)l 
Vi \ dx )\ a sin pi J 


fsinfti a ,,sm 4 pidSi 


A 1 f sin fli u 7 c \ sin 4 7 57i 

2 cos ft t cos fti) |_ y w + ^i) 

J_ ( dV '\ cos (0t + ffi)~ | 

Vi \dx/\ u sm.Pi J 

Pi== Ti (^tXS^pr l )~ T D 


2 sin 3 pi bSj 

.yR cos (0i+^i) bn i 


1 (dVi\ cos (0i -<30 i 

Vi v £*3 AijCOS 2 ^! COS (0i + ftl)J 

p=j_(§Y±\ ( t iw r 2sin8 ^ 

2 Vi \ dx /x u \cos 2 fti A Ly-R COS (<t>i—P\) dn, 

J_/7Vi\ cos (0! + <30 “ j 

Vi \ dx J\ u cos 2 fti cos (0i— fti)J 

T=\A — y ~ 1 - - - • _ 

sin 2 Pi 


dS 

The coefficients Ai and A 3 and the value of -7- ^ have been 

a Mi 

determined for the flow at zero angle of attack and 


dS, 


s lp ~s lpi 


dvi 




sin fl i ~] 


L.cos (0 1 +j8 1 )Jp, 


dSj S ^p i ~ S Jp l 
dni Si~-&i Pi 


(Cl) 


The practical use of equations (59) is identical to the use 
of the corresponding equations (57) for the case of zero 
angle of attack. (See reference 6.) 

CONCLUDING REMARKS. _ 

The method of. characteristics has been applied to bodies 
of revolution at a small angle of attack. Only the first- 
order effects of the angle of attack have been considered. 
The system developed takes into account the effects of the 
entropy variations on the flow phenomena and determines a 
flow that exactly satisfies the boundary conditions within 
the limits of the simplifications assumed. 

The application of the method to practical problems has 
been discussed and two systems are given. The first method 
is numerical and analytical and requires less numerical 
computation but requires the construction of another char- 
acteristic net; whereas the second method is only numerical 
and uses the characteristic net aud some of the numerical 
computations made for the calculations for zero atigle of 
attack. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 22, 1948. 



APPENDIX 


DETERMINATION OF FLOW PROPERTIES AROUND A CIRCULAR CONE AT A SMALL ANGLE OF ATTACK 

Assume a polar coordinate system r,ip,6. Call v T the 
velocity in radial direction, r„ the velocity in normal direction 
to r in the meridian plane <?=Gonstant, and w the component 
normal to the meridian plane (fig. 8) ; that is, 


V r = 


W- 


dr 

dt 

rd\p 

'~w 

r dd sin ^ 


If the phenomenon is conical, 

dv r 

br 

bv n 


dt 


—Q 


br 


=0 


and 


^-0 

br 


52=0 

or 

5i=o 

br 


Therefore, Euler's equations are 

V„ bv r , W bv T V^-\-W 2 


=0 


da* 


b\p r sin ip bd r 
w do* 1 bp , v r Vn— w* cot $ 


r bp 
Vn bw , 


r sin p bd 
bw , 


w 


pr bp 
1 bp 


v r w-r v n w cot p 


r bp r sin p bd rp sin p bd 

The continuity equation is 

2p o r sin p+v* sin^|^+p sin 'P^+ 


(62e) 


r„p cos p-\-w 


bp , bw 
be~^bd z 


(63) 


and the energy equation is 


7 — 1 



Combining equations (62) to (64) results in 


+ 


(62a) 

bw f w~ 

\ 2wVa. 



sin i pbd V a 2 

) a 2 

sin p bd 

=0 (62b) 

wv n f bv n 



*=0 

a 2 \sin pbd 

*p) 



(65) 


For small angles of attack the velocity components can 
be expressed in the form (references 1 and 3) 


V r = Dr 1 +CCVr 2 COS 8 ") 
P B =^ ni +aa* 2 C0S 8 
w=aw 2 sin 8 


( 66 ) 


when the second- and higher-order terms of the angle of 
attack have been neglected. Equation (65) at small angles 
of attack becomes 


,,(2-^)+,.oot #+ ^(i_^) +s ^ 55 _o car 


) 
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The shock is a circular conical shock having its axis inclined 
at an angle 17 with the axis of the cone. The quantities with 
subscript 1 are the quantities corresponding to the case of 
zero angle of attack. (Indeed, the cone is a particular body 
of revolution and, therefore, the considerations made for 
the case of bodies of revolution are still valid.) 

From equations (62a), (62b), and (64b), there results at 
small angles of attack 

1 dp 7 dp = p 

p p b\[/ 


In the meridian plane 0=Constant, therefore, the transfor- 
mation behind the shock is isentropic for small angles of 
attack. If —a&S 2 sin 6 is the variation of entropy in a direc- 
tion normal to the meridian plane 0=Constant (reference 5), 


bS 

sin ^ dfl 


a ASj sin 9 


where A Sa is constant and independent of ^ outside of the 
vortical layer. From equations (626) and (64a) 


— ^AS',=® M ^sin ^+Wa(t? r sin $+v n cos ^)+ 

7 (68) 

If the variation of entropy is small and the term A S 2 can be 
neglected, equation (68) becomes 

w s sin ip= — v rt (69) 

Equation (67) can be written in the following form: 


cot *- 


aWj cos 9 
sin \p 


(70) 


By use of equation (67) and by considering the conditions 
for zero angle of attack, equation (70) becomes 

1 _ 9 


=—v 


cot 'f 


Oil -I 


„ A 4.7-1 cot * A w 2 

Vr 2 I I "r _ 4 9 I / 

2 ^ ^ __ J 9m ^ 


(71) 


Equations (68), (70), and (71) permit the determination 
of the flow around the cone at an angle of attack by means of 
a step-by-step calculation when the calculation for a=0 
has been performed. Consider the hodogTaph plane uv, 
and consider the variation of velocity components v T and v n 
in a meridian plane 0=Constant (fig. 9). Assume that for 
a given value of and 9 a the velocity components v„ v n> and 


V 



w are known. Point P* of the hodograph diagram represents 
the velocity vector OP, corresponding to the velocity at 
every point of the space of coordinate in the plane 
#„= Constant; whereas OQ* represents the values of (©,) p m 
and Q«P a represents the values of (o r )p.- Now, the radius of 
curvature R a of the hodograph diagram is along the line 
Q0P0 and has a value given by (reference 1) 



and, therefore, can be determined from equation (70). At 
P 0 the values of o r , and ©„ t are known from the calculation 
for a=0; therefore, a© rj and av^ can be determined from 
equations (66). Equation (71) can be used in place of 
equation (70) in the following way: The vectors OQ* and 
Q S P 0 represent the values of v H and v T% at P 0 ; the vector OP* 
in the hodograph diagram gives the values of T T j; therefore, 



can be obtained from equation (71). 

Now al any point P* the radius I?*, given from equation 
(70), or the radius R i<t , given from equation (71), is known; 
therefore, from the quantities at P* the quantities at P t ol 
coordinate #»=^*+A^ can be obtained by constructing a 
circle of center C* (where C a P*=i? a or f? 2 j through tho point 
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P 0 until the point P» along the line C„P 6 which is a straight 
line from Ci, is inclined by ft a -{-Aft with the u-axis. There- 
fore, 

(v n )j, + A*=(v„){ cos Aii + (R—v r )f sin A17 
(®r)#+A#=(®»)* sin Aij — (R — v r )f cos Aij -\-R+ 

Inasmuch as the values of v r and i\ at P s have been ob- 
tained, the values of p rj and s, 2 can be determined by 
differences from the values for a=0 with the use of equa- 
tions (66). (If equation (71) is used, the values of and 
0*2 are obtained directly.) With equation (68) the value 

of ^- z can be calculated at P„, and the value of w 2 at P s 

can be obtained. Indeed, A S 2 is constant and has been 
determined from the conditions at the shock. In a similar 
way, all the hodograph diagram can be constructed. If 
necessary, for every point P& a second approximation can 
be determined. 

The calculation of all the flow field must start at the shock. 
For the calculations it is convenient to choose a coordinate 
system having the axis of the conical shock as the axis of 
polar coordinates. In this case, the velocity components 
i>, r„, and w behind the shock can still be expressed in the 
form of equations (66). Indeed, from equations (15) and (19), 


COS »(<.•„,+' d# ‘)2 


= S«^+(a— i7)27* 2t C0S 6 

(72a) 

COS 6 (l^+J 


=V ri +(a—ri)Vr h cos 6 

(72b) 

w,= ac sin d(w 2 — r — v.) 

\ sin ft a Vj 


=(a—i})w ia sin 6 

(72c) 


where v n> , v rt , and w t are the components referred to the 
axis of the conical shock, whereas the components v r6 , «a s , 
and Wt are referred to the axis of the body. Indeed, ij ja 
is constant. 

The calculations start at the shock. After determining 
the flow field for zero angle of attack, the angle of the conical 
shock ft, is known and the velocity components and v, L 
with respect to the axis of the shock for every value of ft are 
also known. In order to determine the flow for the case of 
a small angle of attack, the direction of the undisturbed 
velocity must be rotated at a small angle a — 57 with respect to 
the axis of the shock (fig. 10). The value assumed for a — tj 
fixes the value of a for which the calculations are performed. 
(This value of a is not yet known but is obtained as a result 
of the calculation.) 


V 



For a unitary value of a — rj, the component v Ts behind the 
shock can be determined from equation (14b) {v t =Vt) 
and the component w t , from equation (14c) (w'=w,); 
whereas v„ t can be determined from equation (21a) 
{v n =—Vir") and v N ' is given by equation (14a). The value 
of entropy (a— tj)A S 3 cos 6 can also be determined from the 
equation of the shock, for example, from the difference 
between A S and ASi. When v T) , i\ t , ir 2j , and AS a are known 
behind the shock, all the flow field can be obtained by means 
of equations (68) and (70) or (71). The hodograph diagram 
can be constructed, for example, in the plane B=v. 

The axis u has been chosen in the direction of the undis- 
turbed velocity for zero angle of attack that corresponds to 
the axis of the shock for a=0. For a the undisturbed ve- 
locity has been rotated at a — n with respect to the -M-axis 
(fig. 10) ; therefore, the axis of the shock has not been changed. 
The velocity OP„ behind the shock of figure 9 must be decom- 
posed into (1) a component P a Q„ inclined at ft, and correspond- 
ing to v rt if equation (70) is used or to v T2 if equation (71) 
is used and (2) a component 0 0 Q a corresponding to v„ t or 

In this way, the values of v Tl and v„ L that must be used in 
equations (68), (70), and (71) are the values obtained from 
the calculations for zero angle of attack at the same value 
of ft (that is, ft=ft, for P„) . Because the calculations start 
at the shock, the construction of the hodograph diagram 
must be performed in the direction of decreasing values of ft . 
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At the surface of the body for 0=ir, the component v n 
must be zero ; therefore, when the radius of the hodograph 
diagram passes at the origin of coordinates u and v, the cor- 
responding value of 4 / is equal to i/'o+jj where 4'o is the angle 
of the cone (fig. 10). Because 4 / o is known, the value of n 
and, therefore, of a can be determined. 

The components v r and v n in the plane 0=0 or 0 =t do 
not change when, for the axis of reference, the axis of the 
body is assumed; but the corresponding value of 4/ is in- 
creased at jj (fig. 10). The value of vh changes; the value 
of w ib " can be determined from the value of Wt by means of 
equation (72c). 

For practical calculations it is convenient to use nondi- 
mensional coefficients obtained by dividing all the velocity 
components by the limiting velocity F». The expression 
ajVi can be obtained from equation (39d). 

For small values of a, the values of e rj , and i j/a 

are independent of a and, therefore, the flow for every other 
value of a can be obtained from this determination. The 
calculations can be graphical or analytical. 
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